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Abstract 

A new scheme for relating the absolute value for the kinetic rate constant k to 
the thermodynamic constant K p is developed for gases. In this report the forward 
and reverse rate constants are individually related to the thermodynamic data. 

The kinetic rate constants computed from thermodynamics compare well with the 
current kinetic rate constants. This method is self consistent and does not have 
extensive rules. It is first demonstrated and calibrated by computing the HBr 
reaction from H 2 and Br 2 . The method then is used on other reactions. 

Symbols 

A proportionality constant in equation (16) 

A, classical preexponential factor in forward rate constant 

q species concentration for species i, moles/cc 

□' universal kinetic/thermodynamic link constant -3x10® sec' 1 

E activation energy, calories/mole 

^ temperature constant, equation (16) 

G free energy 

k chemical rate constant 

kk chemical rate constant from classical kinetics 

kt chemical rate constant from thermodynamics 

K,. equilibrium constant based on concentrations 

K equilibrium constant based on equation (8) or pressure 

n 5 P number of moles of i per gram of mixture 

n T total moles of mixture/gram of mixture or 1 /molecular weight 

N constant in modified Arrhenius rate expression eq. (4) 

pj partial pressure of species i, n/ryP atm 

r f forward chemical reaction rate, moles/cc-sec 

R universal gas constant 1.987 cal/mole K or 82.05 atm-cc/mole K 

Sj species i chemical symbol 

t time 

T temperature 

x reaction progress variable, moles reacted per gram of mixture 

X Pease’s concentration, moles/22.4 liters, ref. 6 

A difference operator, products minus reactants 

Copyright © 1996 by the American Institute of Aeronautics and Astronautics, Inc. No copyright is 
asserted in the United States under Title 17, U.S. Code. The U.S. Government has a royalty-free 
license to exercise all rights under the copyright claimed herein for government purposes. All other 
rights are reserved by the copyright owner. 
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p°i standard chemical potential for species i 

v. stoichiometric molar coefficient for species i 

subscripts: 

c value expressed in terms of concentration 

f forward direction 

i for each species i 

p value expressed in terms of pressure 

r reverse direction 

Introduction 

There are several reasons why it would be desirable to relate the kinetic rate 
constants to thermodynamic data. First, modern kinetic models contain many 
reactions with many constants to form large mechanisms.These systems are not 
universal and are complex to assemble from the experimental data base. 
Thermodynamic equilibrium calculations, on the other hand, have been 
simplified and systematized and are only species dependent. Thus, relating the 
kinetic constants to thermodynamics might reduce the work involved in doing 
kinetic calculations. 

Secondly many kinetic constants are still unknown: these could be computed by 
using thermodynamic data; thus, the kinetics calculation could be done without 
performing more experiments. 

The derived relationship that links the chemical rate constants to 
thermodynamics is the ratio equation (1), where the thermodynamic equilibrium 
constant K,. has been established as the forward rate constant k, over the reverse 
rate constant k,: 



Normally one thinks of kinetics as a path process, but the kinetic rate 
constants are state variables like K c . The kinetic constants k are functions 
only of temperature and are therefore state variables and not path variables. 
Equation (1) holds not only at equilibrium, but also at every point of the 
reaction from initial conditions to final state; furthermore, it is independent 
of concentration. 
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The first hurdle to developing an overall relationship between kinetics and 
thermodynamics is to overcome the idea that this is impossible. All of the 
kinetic text books principally state that thermodynamics is only for 
equilibrium calculations whereas kinetics describes the rate of approach to 
equilibrium; thus, they state that thermodynamics would not be useful for 
kinetics calculations. However, as stated, relationship (1) holds for all time, 
not only at equilibrium. Thus, with thermodynamics, the complete energy surface 
is available for ideal gas mixtures for all compositions, not only for 
equilibrium states (from initial reactants to final products). This 
information can be used to compute the kinetic rate constants. The change of 
system free energy with reaction dG/dx is known, and is used to compute the 
kinetic rate constants. 


A second hurdle is the idea that reacting molecules have to go through an 
activated complex energy barrier to complete the reaction. This state may be 
estimated from transition state theory to determine the kinetic constants but 
would be difficult to compute from thermodynamics. But thermodynamics is 
involved in equation (1). Are the activated complex states the same for both 
reaction directions? Does the activated complex state would depend on the 
starting point? 

A connection between thermodynamics and the chemical kinetics constants is 
presented with comparisons to experimental data. 

Description of the Method 

The classical kinetic method can be described for a particular rate by: 

V 1 S 1 +V 2 S 2 =V 3 S 3 +V 4 S 4 > (2) 

where the v f are the stoichiometric coefficients for species Sj. The forward 
rate r f is given by: 

r f = k f C;*C 2 V2 . (3) 

and the kinetic rate constant k, is defined as: 

k f = A f T n exp(-E/RT) (4) 

In the thermodynamic derivation, the kinetic rate is assumed to be related to the 
gradient of the free energy or. 
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, -1 0Gv 

r, a exp( ) . 

f RT5x 


( 5 ) 


The free energy G , ref. 1, is given by: 


G=E(p> i )+£(n i RTln(p i )). 


(6) 


Differentiating G with respect to the reaction progress variable x, and using 
the facts that (dn J /dx=-Av i ); for reactants, (-vjd/dn^ d/dx) ; and for 
products, (y t d/dn { = d/dx) with equation (2) as the model, this results in: 


1 dG 


_ T , =£ -ibr'£ -^ + £ ( v n p fi )-£ ( v n to p ri > . 
RT9x i RT i RT i i 

At equilibrium, the free energy is a minimum, dG/dx= 0, with the equilibrium 
constant Kp given by: 


v n h; 


K = expCE^-E ^)^ 3 ^ 4 


RT 


RT 


„vl v2 

Pi P2 


at equilibrium 


The LHS is the thermodynamic expression for Kp, expressed as: 

„ ,-AG\ 

K P " exp( -if >> 

and is independent of composition or degree of reaction. In addition K,. is given 
by: 


K c =K p (RT)" AVi =K D (RT) (SVfi_SVri) . 


( 7 ) 


( 8 ) 


( 9 ) 


( 10 ) 


At this point, it is important to distinguish chemical kinetic constants, kp and 
k,.. When the kinetic rate is given in terms of concentrations, k,. is used; and 
when the rate is given in terms of pressure p, kp is utilized. 


ii=k c C,V ; p,V ; C ,-A . 

dt c 2 dt p 1 2 1 RT 

k c =k p (RT ) (Ev ”' 1) ; 


Then: 


(ii) 



similar to Kp and Kg in equation (10). 

Returning to equations (5) and (7) and equating the pressure terms (ie. the 
concentration terms) to the concentration terms in equation (3), one obtains: 


k cf =DK p (RT) Sv “ 1 for K < 1. 

where D is a constant to be determined. Note that equation (7) is written for 
the net rate, whereas equation (3) is for the forward rate. Equation (12) 
provides the connection between thermodynamics and the forward chemical kinetic 
constant that shall be used. 

Equation (12) without the inequality cannot hold for both the forward and 
reverse rates because equation (1) would not be satisfied. Remember that 
Kp r =1/Kp f , so when Kp f <1, then Kp r >1. To satisfy the kinetic constant ratio 
equation (1), the reverse rate must be given by: 


k cr =D(RTf'-' for Kp> 1. 


D is an unknown, assumed universal constant. It was determined from the fit of 
the HBr data presented in the next section. It was found to be 3x1 0 8 sec 1 which is 
used throughout this report. 

The results from the thermodynamic method using equations (12) and (13) to 
compute the kinetics for all reactions are compared with experimental data. If 
the chemical potential of the products is greater than the reactants, then 
Kp < 1, and equation (12) is used. If the chemical potential of the products is 
less than the reactants, then the reaction is a cascade reaction with Kp > 1, and 
equation (13) is used. In this formulation only a single constant D and the 
thermodynamic base is sufficient to determine the whole kinetic data base. 

In the next section the hydrogen bromine reactions are discussed, followed by a 
subset of the H/O reactions. In the final section, a random selection of rates 
chosen from reference 4, NSRDS, are presented for comparison. 
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Calibration and Demonstration of Thermodynamic Kinetics 
Method using the Hydrogen-Bromine Reactions 

The HBr system was chosen for several reasons. It was simple, but yet had 
interesting complexities. First there are only 5 species involved and they are 
all in the thermodynamic data set. They are arranged with HBr in the middle as 
follows: 


H Br 

(V /i 


HBr 


1/ \l 

H2 Br2 


Secondly, the three classical rates to describe the system are: 


2 Br 

= Br2 * 

Cl 

A=1.67 E+14, 

n=0., E=-1844. 

Br+I^ 

= HBr+H 

C2 

A=3.4 E+12, 

n=0.5, E=18300. 

HBr+Br 

= H+Br 2 

C3 

A=8.13 E+12, 

n=0.5, E=44500. 


*At low temperatures M is not included, data is from ref. 3. 

These are called “elementary reaction rates,” because kineticists believe they 
represent the specie combinations by which the reaction proceeds. The rates for 
the three reactions are well known and the computed data can be treated like 
experimental data to match with other prediction methods. 

Finally, the overall classical rate for HBr formation is given by: 

d(HBr) = k [HJ[Br 2 p 
dt 5 + [HBr] 

[BrJ 

The rate is inhibited by the HBr product. This is one of the first systems where 
a mechanism could be written down and solved. 

Thermodynamic approach 

In order to obtain the constant D, only one point is required. The formation of 
HBr at 574.3 K , table 2, was used to determine its value at 3 x 10 8 . Individual 
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rate comparisons are presented, followed by the predictions of the 
concentrations. 

The classical mechanism rates above (equation (14)) with the thermodynamic 
predictions could have been used to make concentration predictions. However, a 
generalized equation set was chosen. All possible reactions are written down to 
maintain material balances. In computing thermodynamic equilibrium, all 
species are considered. 

For the thermodynamic method, all six interaction equations are listed below in 
Table 1: 


Table 1. Reactions for t 

he Thermodynamic Hydrogen-Bromine Mechanism 

Equation 

Name 

K p at 574.3 K 

krf co-mol-sec units 
at 574.3 K 

Br 2 =2 Br 

T1 

9.32E-13 

2.90E-4 

H 2 =2H 

T2 

4.10E-35 

1.24E-26 

0.5 H 2 +Br=HBr 

T3 

1.81E11 

6.51E10 

H+0.5 Br 2 =HBr 

T4 

2.73E22 

6.51E10 

H+Br=HBr 

T5 

2.82E28 

1.41E13 

0.5H 2 +0.5Br 2 =HBr 

T6 

1.74E5 

3.0 E8 


One way to compare the above two methods is to use an empirical least squares fit 
of Kp as a function of T, as in the last section of reference 4, in the form: 

Kp = A exp(E;/RT). (16) 

For bimolecular reactions where Kp < 1, k, would equal A(DRT)exp(^/RT) - a form 
similar to the Arrhenius expression, equation (4). Comparing ^ to E gave a hint 
of the desired correlation. This comparison has slightly different numbers, but 
it looked reasonable. * 

A better method is to plot the two relationships versus temperature. A 
comparison of the C2 and RC3 kinetics constants by the two methods are shown in 
figure 1. Equation (12) works well for C2 while equation (13) works for RC3. For 
RC3 the methods agree at lower temperature, but diverge at higher temperature. 

The classical relation may be extrapolated out of its range. 
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Figure 1. Comparison of classical kinetic equations and thermodynamic 
predictions for equations C2 and RC3. 


The two methods for atom/atom recombinations are shown in figure 2. Differences 
between the two methods range from 10' 3 to 10 6 , but the curves follow the same 



Figure 2. Comparison of classical kinetic equations and thermodynamic 
predictions for atom reactions using equation (12) and (13). 
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A better comparison is obtained by plotting the ratio of the thermodynamic HBr 
rates, kt, to classical kinetic rates, kk, versus temperature in figure 3. A 
value of one indicates perfect agreement. Using D of 3.0x10 s , the thermodynamic 
rates are given by equations (12) and (13). 



— * > “ T1:BR2 = 2 BR 
— - C2:BR+H2-HBR+H 

— - C3:HBR + BR = H + BR2 

- .x- - T2:H2 = 2 H 

- • * C1:2 BR = BR2 

— - “ T5:H + BR = HBR 


Figure 3: Ratio of thermodynamic rate to kinetic rate versus temperature 


Rate C2 agrees very well over a wide temperature range. For rate C3, the ratio is 
1/10 the classical; and for the other rates, the ratio at 500 K is 1/100. 500 K is 
the temperature range where the experimental data is compared. The maximum in 
the Br curves is caused by Kp going through one at approximately 1500 K. 

Rather than compare individual equations, it is more intuitive to compare the 
predicted concentrations from the experimental data. The GLSENS program, ref. 

5, is used as the kinetic integrator to compute the concentrations. It can handle 
global reactions and fractional stoichiometric coefficients. 

The conditions of Table 2 were used to obtain the constant D. A value of 3x10® 
resulted in a lower prediction 0.1481 at 870 sec, compared to data of 0.1559, but 
later the predicted results became about 3.5 percent to high. This value 
produced close agreement with the experimental data. The classical method was 5 
percent too low. 
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Table 2: HBr Mole Fraction Data at 574.3 K ; p=1.8 atm; Table XIII of ref. 6. 
Initial Conditions: Xho= 0.657 X Rr9 =0.343 


Time 

seconds 

Data 

^HBr 

Classical 

^HBr 

Thermodynamic 

^HBr 

870 

0.1559 

0.1404 

0.1481 

1470 

0.2295 

0.2165 

0.2312 

2070 

0.2940 

0.2797 

0.3014 

2970 

0.3830 

0.3559 

0.3869 

4770 

0.4877 

0.4613 

0.5052 

5970 

0.5373 

0.5086 

0.5568 

7470 

0.5829 

0.5512 

0.6010 

8970 

0.6102 

0.5814 

- 

Equilibrium 

0.686 




This Table shows a good comparison between both methods and the data. 

Reaction T6 (0.5 H 2 +0.5 Br 2 = HBr is deliberately left out of the reaction set. 

If reaction T6 (molecule/molecule) is included in the reaction set, the 
formation of HBr is very fast with equilibrium being achieved in less than 1 
microsecond. Even though k is small, the concentrations of H 2 and Br 2 are high, 
making the overall rate competitive. Thus, rather than reduce D to include this 
reaction, T6 is left out of the set in order to integrate the rates more like the 
classical set. When T6 was not included, the net conversion agreed with the 
experimental data. The results are shown in Table 2. 

In order to determine the temperature sensitivity of the reaction sets, 

Reference 6 tabulates the experimental k of equation (14) versus temperature. 
The results are tabulated in Table 3. The comparison is good with additional 
roll off with T for the thermodynamic set as T decreases. 

The decrease in HBr formation rate with increasing HBr initial concentration was 
investigated next, using Table XII of reference 6 for data. The effect of 
increased initial HBr concentration is shown in Table 4. The experimental data, 
given in units of moles/22.4 liters of H2 was converted to moles/cc of HBr by: 

9 

CHBr( m °l es / cc ) =X(moles H^/22.- 4 liters)* 

22,400cc 





































Table 3: Predicted Temperature Sensitivity of Reaction Sets 
and Comparison with Table XIV of ref. 6. 

Initial Conditions : Xh 2 =0.657, X^O.343, time=14.5 min 


Temperature K 

D 

'atm 

^experimental 

^classical 

Khermodynamic 

574.3 

1.806 

25.62 

23.06 

24.56 

550.5 

1.731 

4.760 

4.862 

4.186 

524.4 

1.649 

0.7782 

0.7386 

0.5058 

498.7 i 

1.568 | 

0.0970 

0.0950 1 

0.0510 


* values of k are (cc/mole) 0 5 min 1 


Table 4: Effect of Increased Initial HBr Concentration at 574.3 K and 25 min. 

Initial H2 and Br2 = 2.009E-05 


Initial 

HBr* 

Patm 

Gain in 
Cone. HBr 

Classical 

HBr T 

Classical 

AHBr 

Thermo 

HBr x 

Thermo 

AHBr 

0 . 

1.89 

1.066E-5 

8.543E-6 

8.543E-6 

9.20E-6 

9.20E-6 

4.1 IE-5 

3.83 

9.540E-6 

4.81 E-5 

7.0E-6 

5.49E-5 

1.38E-5 

5.58E-5 

4.53 

8.18E-6 

6.25E-5 

6.7E-6 

6.97E-5 

1.39E-5 

7.55E-5 

5.46 

7.20E-6 

8.17E-5 

6.2E-6 

8.94E-5 

1.39E-5 

8.62E-5 

5.96 

8.25E-6 

9.21 E-5 

5.9E-6 

1.00E-4 

1.38E-5 


+ all concentrations are given in moles/cc 


Note that as HBr is added to the mixture, the pressure increases. The classical 
equations overpredict the suppression of HBr while the thermodynamic method 
predicts an increase, but is not sensitive to the quantity of HBr added. The 
dilution of the reactants with HBr is compensated for by the increase in 
pressure, keeping the concentrations of H 2 and Br 2 fixed. The predicted final 
concentration of HBr is within 10 percent of the data for both methods. 
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The final comparison examines the effect of varying the initial H 2 and Br 2 
concentrations, using the data from Table XI of ref. 6. The results are shown in 
Table 5. 


Table 5: Effect of Changing Initial Concentrations of H 2 and Br 2 , 
Temperature 574,3 K 


Initial 

h 2 

Mole 

Fraction 

Initial 

Br 2 

Mole 

Fraction 

Pressure 

atm 

Time 

minutes 

Exper 

HBr 

moles/cc 

Classic 

HBr 

moles/cc 

Thermo 

HBr 

moles/cc 

0.657 

0.343 

1.805 

14.5 

5.97E-6 

5.38E-6 

5.65E-6 

0.655 

0.345 

0.9252 

19.5 

2.87E-6 

2.67E-6 

3.78E-6 

0.380 

0.620 

1.718 

15.0 

4.39E-6 

4.13E-6 

1.01 E-5^ 

0.378 

0.622 

0.8638 

15.0 

1.70E-6 

1.52E-6 

5.12E-6^ 

0.5 

0.5 

1.010 

20.0 

3.14E-6 

2.9E-6 

6.11E-6^ 

0.5 

0.5 

3.787 

25.0 

2.44E-5 

2.24E-5 

2.75E-5 

0.5 

0.5 

2.840 

15.0 

1.10E-5 

9.95E-6 

1.34E-5 

0.5 

0.5 

1.893 

15.0 

5.98E-6 

5.56E-6 

8.94E-6 


★ Indicates 


large deviation from experimental data 


Both the classical and thermodynamic predictions agree well with the 
experimental data, except for the three starred cases where the H 2 was low. For 
the thermodynamic kinetic setT3 (0.5 H 2 + Br=HBr) is in equilibrium and is 
waiting for T 1 (Br 2 =2 Br) for bromine. The time constant for T1 at the low 
temperature of 574.3 K is 60 minutes (1./k f ). So HBr can be expressed as: 


HBr = H 2 1/2 Br ; 


dBr _ dBr, 

dt dt 


dHBr 

dt 


The major difference in deriving mechanisms is that the classical method assumes 
that Cl (2Br=Br 2 ) is equilibrated, whereas in actuality it is the rate 
controlling step and the rate of formation of HBr (T3) is equilibrated. 


This comparison shows that the thermodynamic method is promising for predicting 
reacting performance. However more examples are required to show that the 
results are not a peculiarity of the HBr system. 





























































Demonstration of Thermodynamic Kinetics Method 
on the Hydrogen/Oxygen System. 

To determine how the thermodynamic method worked on another system, the first 
six rates of the H/O system from a recent publication, reference 7 are compared. 
The ratio of thermodynamic to classical rates are shown in figure 4. 



Figure 4. Ratio of thermodynamic kinetics constant to classical kinetics data 
for some reaction of the Hydrogen/Oxygen system. 


At higher temperatures, disagreements were on the order of one order of 
magnitude; however, at the low temperature of 500 K, the discrepancy is large - 
10 4 for the 0+H 2 reaction. This is because the free radical concentration is low 
at low temperatures, making measurement accuracy low. Also at low temperature 
the rates are low so that systematic errors are magnified on a ratio scale. The 
recommended temperature range for this data is 1400 to 1900 K. The accuracy of 
the data is estimated to be 13 percent, whereas the deviation of the 
thermodynamic prediction is an order of magnitude, at high temperature. 

The error decreases as the temperature increases, approaching a factor of 10 at 
2500 K. Again below 1000 K, the classical kinetic data is an extrapolation of 
high temperature which may account for the large deviation and errors look 
bigger. These thermodynamic predictions were made with no change in D (3x1 0 8 ). 
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Demonstration of the Thermodynamic Kinetic Method 
on a Random Selection of Rates 

As a final comparison, the thermodynamic kinetic method is compared to the 
classical method for a random selection of kinetic rates. Rates from reference 
4, five with a thermodynamic constant greater than one, Kp> 1, and then five 
reactions where Kp < 1 were selected. The first rates listed were with Kp > 1, 
which turned out to involve the hydrogen atom H. These are shown in figure 5. 



Figure 5. Ratio of thermodynamic rate to kinetic rate for Kp > 1 , ref. 4. 

When Kp > 1, the rates are given by equation (13) or for a bimolecular reaction 
the rate is k* = DRT. When Kp > 1, the rate constant is independent of the initial 
chemical potential, as long as the final chemical potential is below the initial 
value. The agreement is very good, especially for the ozone 0 3 and the peroxide 
H0 2 reactions. 




The comparison with Kp < 1 is shown in figure 6. 



Figure 6: Ratio of thermodynamic rate to kinetic rate for Kp < 1. 

Comparing the deviations in figures 5 and figure 6 shows that about the same 
trends occur with large errors at low temperatures and decreasing as the 
temperature increased. The value of D is still 3x10®. 

In general the thermodynamic predictions show the same trend as the classical 
kinetic rates. More comparisons need to be carried out to extensively validate 
the thermodynamic method for kinetic rate prediction. 

Conclusions 

Using thermodynamic data to compute chemical kinetics constants may be a viable 
alternative for gas phase reactions. The value of the proportionality constant 
D appears to be universal and equal to 3x10® for the cases studied. Changing the 
value of D would only shift the predicted kinetic rate constant higher or lower 
for both the forward and reverse rate for the whole temperature range. 


The thermodynamic method is given by: 




(12)&(13) 


k c =DK p (RT f vl for K p < 1 or 

k c =D(RT) SVfi 1 for K p > 1 

where D=3xl0 8 sec -1 . 

No specific mechanism is used. All of the species have to interact according to 
their standard free energy level, p°. All species have a chance to interact like 
in the equilibrium code. 

There is a strong corellation between the two methods. It is remarkable that for 
the reactions studied, this simple thermodynamics kinetics model works so well 
and is consistant for a variety of systems. 
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